Bi-Polaron and iV-Polaron Binding Energie^] 
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The binding of polarons, or its absence, is an old and subtle topic. Here we prove two things 
rigorously. First, the transition from many-body collapse to the existence of a thermodynamic limit 
for N polarons occurs precisely at U = 2a, where U is the electronic Coulomb repulsion and a is 
the polaron coupling constant. Second, if U is large enough, there is no multi-polaron binding of 
any kind. Considering the known fact that there is binding for some U > 2a, these conclusions are 
not obvious and their proof has been an open problem for some time. 

PACS numbers: 03.65.-w, 71.38.-k, ll.10.-z 



The (large) polaron, first considered by H. Frohlich [l[ 
in 1937, is a model of an electron moving in three dimen- 
sions and interacting with the quantized optical modes 
of a polar crystal. In suitable units, its Hamiltonian is 



=p 2 



o t (k)o(k) dk 
1 



(1) 



V2n 



[a(k) exp(ik • x) + h.c] dk , 



where a(k) are the annihilation operators of the scalar, 
longitudinal modes (with [a(k), a^(k')] = S(k — k')), p 
is the momentum of an electron, and a is the coupling 
constant. (Other authors have used a different conven- 
tion, where a is replaced by a/^/2 [H, In the ground 
state, with energy E^\ the electron accompanied by the 
localized excitations of the phonon field constitutes the 
polaron. Through the years, the polaron has served both 
as a model for an electron in an ionic crystal and as 
a simple model for a dressed particle in nonrelativistic 
quantum field theory. 

Of great physical interest is the binding energy of N 
polarons, with Hamiltonian 



H u )= Y,P 2 j+ I « f (k)a(k)dk 
fa N f 1 



(2) 



u Yl i^-^r 1 

l<i<j<N 



and ground state energy E- 



(N) 



Here, U > is the 



Coulomb repulsion parameter, equal to e . The deriva- 
tion of in [1] implies that U > 2a, and this is crucial 
for thermodynamic stability, as we shall see. 

We first consider the bipolaron binding energy AEjj = 
2E^ — E^p . For some time this was thought to be zero 
for all U > 2a, on the basis of an inadequate variational 
calculation, but it is now known Q to be positive for some 



U > 2a. The question we address is whether AEjj — 
for U sufficiently large. We will show that there is a fi- 
nite constant C such that AEjj = if U/a > C. It is 
understood that the effective interaction induced by the 
phonon field for two polarons at large distances d is ap- 
proximately Coulomb- like —2a/d, but this alone does not 
preclude binding. The known existence of bipolarons for 
some U > 2a is an effect of correlations. It is a priori con- 
ceivable that correlations lead to an effective attraction 
that is stronger than Coulomb at large distances. If it 
were, for example, equal to {2a/d) log(log(log(d))), then 
this minuscule perturbation of Coulomb's law, which 
would be virtually undetectable by a variational calcu- 
lation, would result in binding for all U. The finiteness 
of C is a problem that has resisted a definitive resolution 
for many years. 

The second problem we consider is the existence of the 
thermodynamic limit. For large N, physical intuition 
suggests that E^' ~ —const N. This supposition is 
known to be false if U < 2a. Indeed, it was shown in Q 



(AO 



that, even with the Pauli principle, E [ 
when U < 2a. Absent the Pauli principle, E\? '' would 



-const N 7 / 3 

(AO 



behave even worse, as —const iV 3 . It is also known [§| 
that E { y ] > -const N 2 if U > 2a. The latter bound 
ought to be —const N instead, and we prove this for all 
U > 2a. Even more is true; there is a number U c (a) such 
that when U > U c (a), then E^ = NE^\ i.e., there is 
no binding whatsoever. There will, of course, be an in- 
termediate region in which bound complexes form, a gas 
of bipolarons, for example, or a crystal. The Pauli prin- 
ciple plays no role in our considerations, and our results 
hold equally for fermions and bosons. 

The following rigorous results concerning EW will be 
important in our analysis, (i) For all a, E^ < —a 
[J]. For small a, E^ ~ —a according to the lower 
bound in [Hj|, which is E^ > —a — a 2 /3. (ii) For all 
a, E^ < -C P a 2 0, where C P = 0.109 is the number 
determined by Pekar's integral equation 0. (iii) As^ 
totically, as a — > oo, E^ ~ —Cpa 2 according to [a, |9J]. 
(iv) There is a representation for Eq N ^ in terms of path 
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integrals. In terms of the partition function Z^ N \T) 
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limr^ooT'- 1 logZW(T). 



(Strictly speaking, Z^ N \T) does not exist because of the 



translation invariance of H, 



(A f ) 



and the infinite number 



of phonon modes. These technicalities can be handled by 
inserting appropriate cutoffs, to be removed at the end of 
the calculation 1(J Tlj].) ^ was shown in 12| that after 
one integrates out the phonon variables, Z^ N \T) has a 
functional integral representation 



AN) 



(T) = / d//W exp 
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Jo 



e-^dtds 
|Xj(t) -Xj-(s)| 



(3) 

where d^S N ^ is Wiener measure on all T-periodic paths 
(xi(t), . . . , Xjv(t)). (Strictly speaking, i — s has to be 
understood modulo T, but this is irrelevant as T — > oo.) 

We shall now state our results as three theorems, and 
sketch their proofs. While we postpone the discussion of 
technical details to a subsequent paper, the full structure 
and concept of the proofs are visible in this letter. We 
use the symbols ci,C2,... to denote various calculable 
positive constants that arise in the proof. Our results 
are strongest in the bipolaron case, which is also easier 
than the general case and illustrates our concepts most 
clearly. 

Theorem 1 (Absence of bipolaron binding). 

There is a computable constant C, such that if U > Ca, 

(2) 

then the energy expectation [$> , Hjj ' for any bipolaron 
wave function \f is strictly bigger than 2E^ 1 \ i.e., there 
is no binding if U > Ca. 

It is convenient to structure the proof in four steps. 

Step 1. Partition of the interparticle distance: We fix 
a length I, whose value will later be chosen proportional 
to a^ 1 , and partition the relative distance r = |xi — X2I 
between the particles into spherical shell-like regions of 
radial size 2 k ~ 1 £ < r < 2 k £ with k = 1,2,.... This 
partitioning is one of the key points of our analysis. In 
addition there is the k = region, where the particle 
separation is between zero and £. Because of the uncer- 
tainty principle these regions have to overlap a bit, but 
this can be easily handled and we ignore it for for sim- 
plicity. There is a kinetic energy cost for localizing the 
particles according to this partition, which is a2~ 2k £~ 2 
in the shell k. In the next step we look at the energy of 
the particles localized to one of these shell-like regions. 

Step 2. Further localization for well- separated parti- 
cles: For k > 1 we further localize the particles into indi- 
vidual boxes of size 2 k ~ 3 £. This costs another localization 
error C22~ 2k £~ 2 . Because the separation exceeds 2 k ~ 1 £, 
the two particles cannot be in the same or neighboring 
boxes. From the path integral ([3]), but now with the 
Xj(t)'s constrained to their respective boxes, we see that 
the separated particles feel an effective Coulomb-like at- 
tractive potential. However, this can contribute at worst 
— c 3 a2~ k £~ x to the energy. But the Coulomb repulsion 



is at least U2 k £ 1 , which implies that the total energy 
exceeds 2E^ if 



U2- k r 1 > c 3 a2' k £r 1 + (ci + c 2 )2- 2k t 



(4) 



If this inequality holds for k = 1, it holds for all k > 2 as 
well. Thus, if we can deal with the k = region, we will 
establish that binding is not possible if 



Ua- 1 > c 3 + (ci + c 2 )/{2£a) . 



(5) 



Step 3. The region of no minimal separation: In the 
k = region, the Coulomb repulsion is at least U^ 1 , 
but, since there is no minimal separation, we have no 
direct handle on the possible attraction due to the field. 

We need a lemma, which we will prove in Step 4. It con- 

(2) 

cerns Eq , the energy of the bipolaron with no Coulomb 
repulsion, i.e., U = 0; 

4 2) > 2E^ 



7a 2 /3 for all a . 



(6) 



Assuming this, the total energy in the k = region ex- 
ceeds 2E^ provided 



Ur 1 > 7a 2 /3 + Cl £- 2 , 
that is, no binding occurs if 

Ua- 1 > 7£a/3 + Cl /{£a) . 



(7) 



(8) 



Setting the right sides of ([5]) and ((8|) equal leads to the 
choice £ = c^/a and to absence of binding if U > Ca, as 
asserted. 

Step 4- The universal lower bound ([6]).' In this step, 
U = and we denote the a-dependence of energies ex- 
plicitly. We first note that 



E^{2a) > 2E i - 1 \a) - 4a 2 /3 . 



(9) 



This follows from the lower bound E^(a) > —a — a 2 /3 
in [5( and the upper bound E^\a) < —a in [4|, stated 
above. So ((6|) will follow if we can prove that 



E { 2 \a) > £«(2a) 



(10) 



For this purpose we go back to the functional integral ([3j> 
and use Schwarz's inequality (e a+fc ) < (e 2a ) 1 / 2 (e 2fc ) 1 / 2 , 
where (•) now denotes expectation with respect to Wiener 
measure. We choose a to be the sum of the two terms 
i = j = 1 and i = j = 2 in ([3]), and b to be the 
mixed terms i ^ j. Since (e 20 ) 1 / 2 ~ e _TB< )< - 2a - 1 for large 
T, inequality (|10p will be achieved if we can show that 
(e 26 ) 1 / 2 ~ e Ta . At first sight, the double path integral 
(e 2b ) looks like that for a positronium-like atom, i.e., two 
particles attracting each other through a Coulomb force 
with coupling constant 4a. The trouble is that the inter- 
action in is at different times, i.e., \xi(t) — ^(s) A 



3 



simple application of Jensen's inequality, however, shows 
that we can fix the time difference u = t — s and bound 



(e 2b )< 



dp} 'exp 



4a 



dt 



o |xi(i)-x 2 (t-«)| 



Because of the T-periodic time translation invariance of 
the Wiener measure, the path integral is, in fact, inde- 
pendent of u. Hence we get the positronium-like answer 
as a bound. This completes our argument for the univer- 
sal bound ©, and hence the absence of bipolaron binding 
for sufficiently large Ufa. 

We now return to the case of general particle number 
N. As we noted above, there is no thermodynamic limit 
if U < 2a, even with Fermi statistics Q. In contrast, 

Theorem 2 (Thermodynamic stability for the 
iV-polaron system) . If U > 2a the energy of N parti- 
cles (polarons) is bounded below as 

E { u ] > -C(U,a) N 

where C(U,a) is finite and independent of N, but can 
depend on U, a and on statistics. 

Our upper bound on C(U, a) goes to +oo as U goes 
down to 2a, but we are not certain that this divergence 
reflects the true situation. 

Theorem 2 implies the existence of the thermodynamic 



limit, lim 



iV-VOO 



-1E?W - 



N- X E, 



u 



C'{U,a), when U > 2a. 
The reason is that is sub-additive, i.e., e\^ +M ^ < 

E^ +i?y M ' ) (by considering N particles in a ground state 
located near Princeton and M particles in a ground state 
located near Charlottesville). This fact, together with 
the linear lower bound from Theorem 2, implies the ex- 
istence of the thermodynamic limit, see [l3], Sec. 14.2]. 

Theorem 2 is an essential ingredient for our proof of 
Theorem 3 about the absence of any binding for large 
enough U. 

The proof of Theorem 2 does not use the partition- 
ing and localization of Theorem I. Instead we bound 
from below by the average over translations of a 
Hamiltonian pertaining to a finite-size box and with a 
short range, i.e., Yukawa-like interaction. This 'sliding 
method' of localization was introduced in [3] , and later 
used in [TBI, to analyze Foldy's law for bosonic jellium. 

The localization is accomplished by choosing a func- 
tion x(x) with finite range L and integral J x(x) 2 cix = 
1. Next, consider the function /(x) = (I — e~"l x lx * 
X(x))/|x|, where * means convolution. Essentially, / 
is the difference of the Coulomb potential and a cut-off 
Yukawa potential. If uj is large enough, / will be positive 
definite, i.e., have a positive Fourier transform, as shown 
in [13]. The crucial inequality, then, is 

dx ebd (£. 8(x - x 4 ) - pt (x )) /(x _ x ') 

x (E/( x '- x i)-^ x ')) >° (11) 



for any x-dependent operator p(x). We apply this to 
V2tt 



P(x) 



J ke tk *a(k)dk. 



In the classical jellium case, p = p 1 is just the background 
density, whereas here it is a fluctuating quantum field. If 
we multiply out the various terms in the resulting 
inequality can be written as 



H 



(N) 



> L~ 3 [ dzH z + (U- 2a) V |x 

i<3 



where H z is a Hamiltonian of the particles that happen to 
lie in a box centered at z and with sides of length L. It is 
crucial to note the fact that /(0) is finite, which allows us 
to replace the apparent, unwanted Coulomb self-energy 
in (jlip by a term linear in N. 

The particles interact through Yukawa-like forces. The 
field appears in H z as p(x)x(x — z), so it is localized to 
the box as well. Thus H z refers to a quantum-mechanical 
problem confined to a box of a fixed size but with an 
indeterminate number n z of particles. 

Once the localization to finite-size boxes is established 
we can follow the analysis in [3J , based on the commuta- 
tor bounds in Q , to show that H z is bounded from below 
by — c^a 2 n\ — D{a, L)n z , where D{a, L) — > oo as L — > 0. 
On the other hand, 



Xj| 1 > cgL 4 / dzn z (n z — I 



El* 

i<j 



For any U > 2a, the length L can thus be chosen small 
enough such that (U — 2a)ce/L > C5. The price paid 
for this is the energy —D{a,L)N which can be large, 
but is finite. This concludes the proof of thermodynamic 
stability. 

As U increases from 2a the system is thus stable but 
can form many-body bound complexes such as the bipo- 
laron - perhaps even a periodic super-lattice. This is a 
largely unexplored area. Eventually, no binding is possi- 
ble, as the following theorem asserts. 

Theorem 3 (No binding for large U). There is a 
computable constant U c (a) such that if U > U c {a), then 
the ground state energy equals E { ^ ] = NE™ for all N. 

Although the optimum value of U c (a) might depend 
on particle statistics, our bound does not. We can prove 
that U c (a) < const a for large a. We believe, but we 
cannot prove, that this linear law holds for all a. 

Step 1. We use a similar partitioning as in Theorem 1, 
but relative to nearest neighbor distances. That is, each 
Xj is localized in some shell-like region 2 ki ~~ 1 £ < t» < 2 ki £, 
where ti is the distance between particle i and its nearest 
neighbor in the configuration (xi, x 2 , . . . , xjy). 

Then, as in Step 2 of Theorem 1, we localize each par- 
ticle i in a box whose size is smaller than, but comparable 
to, 2 ki £. As in the bipolaron case we have to remember 



4 



the kinetic energy associated with this two-fold localiza- 
tion. It is a geometric fact that any Xj can be the nearest 
neighbor of at most 12 other particles. This allows us to 
compensate the localization energy by sacrificing part of 
the Coulomb repulsion. 

Step 2. With every particle thus localized in some 
box we write the functional integral for the ground state 
energy as in ([3]) , except that we now include the Coulomb 
repulsion as well as the polaronic attraction terms. The 
exponential now contains 



N 

a x 
2 ^ 



T ,T e -\ts\ dtds 



\Xi(t)-Xj(s) 



l<i<j<N 



dt 



o \Mt)-Xj(t)\ 



We relabel the particles so that ki = ... = km = and 
fc/vf+i, • ■ • , k,N > 1, i.e., such that particles 1,2, . . . , M are 
precisely those having a nearest neighbor within distance 
I. Accordingly, we split the sum in the exponential into 
three pieces. The first one corresponds to the total energy 
of the M particles with ki = 0, 



M 

a x -\ 

9 2.^ 



T ,T e -\t-s\ dtds 



i,j=l J ° J ° 



\Xi(t)-Xj(s) 



- u E 

l<i<j<M 



dt 



o l x i(*) - x i(*)l 



The second piece corresponds to the polaronic self-energy 
for particles i > M, that is, 



a 
2 



N 

E 

i=M+\ 



dtds 



Jo Jo l*i(t) -Xi(»| 



The third is 




t e -\t-s 



dt ds 



-U 



o Jo — Xj(s)| 

T dt S 



o l x i(*) - x i(*)l 



If we keep in mind the confinement of the particles 
to their individual boxes, and how the distances between 
these boxes are related to their sizes, we see that the third 
piece is necessarily negative provided U/a is large. This 
condition is independent of the parameter I. We are left 
with the first and second piece. The second just gives us 
the energy (N — M)E^> after integration. For the first 
piece we write U = U\ + U2 with U\ > 2a. Since the U2- 
part of the Coulomb repulsion is bounded from below by 
const U2M/I by construction, Theorem 2 shows that the 
total energy of the first piece is bounded from below by 
—C(Ui, a)M + const U2M/I. This energy will be bigger 
than the sum of ME$) and the (^-dependent) localiza- 
tion error, provided U2 is large enough. This completes 
the proof of Theorem 3. 



The authors are grateful to Herbert Spohn for making 
us aware of the problem of proving absence of binding for 
large U. Grants from the U.S. National Science Foun- 
dation are gratefully acknowledged: PHY-0652854 (EX. 
and R.F.), PHY-0845292 (R.S.). 



[5] 
[6] 
[7] 
[8] 
[9] 

[10] 
[11] 

[12] 
[13] 
[14] 

[15] 



[1] 

[2] 

[3] 
[4] 



© 2010 by the authors. This paper may be reproduced, 
in its entirety, for non-commercial purposes. 



Electronic address 
Electronic address 
Electronic address 
Electronic address 



: rlfrank@math.princeton.edu 
: lieb@princeton.edu 
: rseiring@princeton.edu 
: let@virginia.edu 
H. Frohlich, Theory of electrical breakdown in ionic crys- 
tals. Proc. R. Soc. Lond. A 160, 230-241 (1937). 
J. T. Devreese, F. M. Peeters, G. Verbist, Large bipo- 
larons in two and three dimensions. Phys. Rev. B 43, 
2712-2720 (1991). 

M. Griesemer, J. Schach M0ller, Bounds on the min- 
imal energy of translation invariant N-polaron sys- 
tems. Comm. Math. Phys., to appear. Preprint (2009), 
larXiv:0907.0752V l. 

T-D. Lee, D. Pines, The motion of slow electrons in polar 
crystals. Phys. Rev. 88 (1952), 960-961; M. Gurari, Self- 
energy of slow electrons in polar materials. Phil. Mag. 
Ser. 7 44:350 (1953), 329-336; T. D. Lee, F. Low, D. 
Pines, The motion of slow electrons in a polar crystal. 
Phys. Rev. 90, 297-302 (1953). 

E. H. Lieb, K. Yamazaki, Ground-state energy and effec- 
tive mass of the polaron. Phys. Rev. Ill, 728-722 (1958). 
S. I. Pekar, Untersuchung iiber die Elektronentheorie der 
Kristalle, Berlin, Akad. Verlag (1954). 
S. J. Miyake, Strong coupling limit of the polaron ground 
state. J. Phys. Soc. Jpn. 38 (1975), 181-182. 
M. Donsker, S. R. S. Varadhan, Asymptotics for the po- 
laron. Comm. Pure Appl. Math. 36, 505-528 (1983). 
E. H. Lieb, L. E. Fhomas, Exact ground state energy 
of the strong-coupling polaron. Comm. Math. Phys. 183, 
no. 3, 511-519 (1997). Erratum: ibid. 188, no. 2, 499-500 
(1997). 

G. Roepstorff, Path integral approach to quantum phy- 
sics. Berlin-Heidelberg-New York, Springer, 1994. 

H. Spohn, The polaron functional integral. In: Stochas- 
tic processes and their applications, Dordrecht-Boston- 
London, Kluwer, 1990. 

R. P. Feynman, Slow electrons in a polar crystal. Phys. 
Rev. 97, 660-665 (1955). 

E. H. Lieb, R. Seiringer, The Stability of Matter in Quan- 
tum Mechanics, Cambridge (2010). 

J. G. Conlon, E. H. Lieb, H-T. Yau, The N 7/5 law for 
charged bosons. Comm. Math. Phys. 116 (1988), no. 3, 
417-448. 

E. H. Lieb, J. P. Solovej, Ground state energy of the one- 
component charged Bose gas. Comm. Math. Phys. 217 
(2001), no. 1, 127-163. Erratum: ibid. 225, no. 1, 219- 
221 (2002). 



